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Abstract 

We consider finite harmonic chain (consisting of A'^ classical particles) 
plus dissipative force acting on one particle (called dissipating particle) 
only. We want to prove that "in the generic case" the energy (per particle) 
for the whole system tends to zero in the large time limit t — >■ oo and then 
in the large A*' limit. "In the generic case" means: for almost all initial 
conditions and for almost any choice of the dissipating particle, in the 
thermodynamic limit N ^ oo. 

1 Introduction 

The energy H{t) in Hamiltonian systems is conserved. After adding dissipation 
terms, H(t) is not conserved anymore but becomes a non-increasing function. 
We consider finite harmonic chain (consisting oi N < oo classical particles) plus 
dissipative force acting on one particle (having number n and called dissipating 
particle) only. It appears that H{N, n, t) is very irregular with respect to N and 
n. However, we want to prove that in the general case the energy per particle 
N~^H{N,n,t) tends to zero in the limit where first t oo and then TV — ;> oo. 
"In the general case" means: for almost all initial conditions and for almost any 
choice of the dissipating particle, see exact definitions below. 

Finite linear hamiltonian systems (including also additional linear terms) 
were studied by many authors, see for example [H 121 IH [S], but with the goals 
different from ours. The papers [SlIZl are closer to ours, in [7] the authors study 
the cases when the dimension of Lq (see next section) is zero. 

Our main formulas concern linear algebra and our main conventions about 
notation are as follows: all vectors are row vectors, there will be introduced two 
scalar products {,)k,k = 1,2, and we denote the linear span of the set R of 
vectors by {R). 

2 Main results 

We consider the system of N particles with the phase space 

L = L2N= = {V' = {q,p) ■■ g - (gi, . . . , qN),P = {Pi, . . . ,pjv) e R^}, 
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with the scalar product 

N 

(V',V')i = 13(ft?i+PiPi) 

and with quadratic Hamiltonian (energy) 

1 ^ 1 

Hii;) = T + U, t=-^pI u=-Y,v{i,j)mo, 

i=l i,j 

where V = {V{i,j)) is a (symmetric) positively definite {N x A^)-matrix. We 
assume unit masses of the particles. The phase space is the direct sum of the 
orthogonal coordinate and momentum subspaccs, Z^-* and Z^'* correspondingly, 
with the induced scalar products {q,q')i and {p,p')i. 

The dynamics of this system is defined by the following system of equations, 
i = l,...,N: 



dt 



dt 

d^ _ _dU 
dt dqi 



- apiSin = - ^ V{i,j)qj - apidin, 



where we appended a dissipation force for one particle ONLY, that is, we fix 

particle with number n and assume that it is subjected to dissipation with fixed 
a > 0. More generally, such system can be written in the matrix form 

ijj = Aij;, ip € L, 

with {2N X 27V)-matrix 

with the matrices V,D and the unit matrix E acting in R^. The dissipation 

matrix D ^ is symmetric and non-negative, in our case it is diagonal with the 
only non-zero element. The dynamics can also be presented as the equivalent 
linear second order ODE system in K^: 



q + Dq + Vq = 0, g = (^i, . • . , gjv) e 



r>N 



For any initial ^{0) G L the solution is tp{t) = e*^^(0). It is well-known (and 
easy to check) that the energy is non-increasing and 

^Him)--iDp,ph. (1) 

Definition 2.1 Lo C L is the subset of elements ip of the phase space L for 
which the energy is conserved, 

Lo = {V e L : ^if(e*^V) = 0, > o}. 
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L_ C L is the subset of elements tp of the phase space L for which the energy 
tends to zero as t ~¥ oo, 

L_ = {i) eL: Hie^'^ij) 0, t ^ oo}. 

Theorem 2.1 Subsets Lq and L^are orthogonal linear subspaces of L, invari- 
ant with respect to the dynamics, and L is their direct sum, 

L = Lq®L^. 

If ip ^ L-j then H^e^^tp) ~^ exponentially fast. More exactly, for some 
constants ci,C2 > (ci — ci{ip) depends on ip), 

H{e*^ip) < ci exp(-C2t) 

Our main result concerns the chain of harmonic oscillators with the Hamil- 
tonian 

1 ^ W ^ Af 

fc=l k=l k=2 

where a;o,wi,a > 0. In terms of the deviations 

qk = Xk - ka 

the Hamiltonian becomes 

N N N 

k=l k=l k=2 

We shall denote by gcd(ai, 02) the greatest common divisor of the natural num- 
bers fli, a2. 

Theorem 2.2 The dimension Dn{N) = dim_Lo.„ enjoys the properties 

Dn{N)=gcd{N,2n^l)-l. (2) 
For any e > there exists some c{e) such that 

N 

N 



1 ^ 

S{N) = -J2D,,{N)<c{e)N^. (3) 



as N ^ 00. 

For some absolute constant c > and any Nq > 0, 



No ^ N 

N^No 
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The following examples show that the function _D„ ( A'^) is very irregular in n 
and N: 

• If = 2'= for some fc > 0, then Dr,{N) = for ah n. 

• If = 2to - 1 and n m, then Dn{N) ^ 2m ~ 2 = N - 1. 

Note that S{N) is the mean dimension if the particle n is chosen randomly, and 
N~^S{N) is the mean energy per particle. Theorem l2.2l shows that if the initial 
energies of the particles in the system are uniformly bounded as — )> oo, and 
thus the initial energy of the whole system is not more than of order A^ (that 
is, proportional to the number of particles), then the energy remaining in the 
system after a long time is very small, namely of the order A^"^, and the energy 
per particle is of order A^^^^*^. It can be said, roughly speaking, that in the 
thermodynamic limit N ^ oo the remaining mean energy per particle is zero 
(even if the kinetic energy is pumping out only of a single particle) . 

3 Proofs 

3.1 Proof of Theorem 

Define the second scalar product in L by 

(V',^')2 = {Vq,q')i + {p,p')i, 

and denote by e„ = (0, . . . , 0, 1, 0, . . . , 0) € ij^^ the coordinate vector with 1 at 
the nth place. Define also the vector g„ = (0, e„) G L. 
We shall prove that 

io - {V- e L : (e*^^, g„)i = 0, > 0} = {V £ L : (e*^V, 9n)2 = 0, > 0}. 

(5) 

Indeed, we have 

This is just rewriting of the right-hand side of ([Ij, 

[Dp,p)i = apl{t), p^{t) = (e*^V,ff«)i- 

Linearity of Lq follows from ([S]) , the invariance of Lq follows by differentiation 
of the equality (6*^-07 = 0. 

Define a subspace ly C as the linear span of the vectors V^Cn, j ~ 
0, 1, . . ., and introduce the subspaces L^'') as the orthogonal complements to Lq 
with respect to the scalar products (, )fe, = 1, 2 correspondingly. 

Lemma 3.1 

. 1. L(i) =i(2). 
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. 2. = e L : q,pelv}. 

• 3. Z/'^^ is invariant with respect to A. 

Proof. We will use the following identity 

A = IQ- aV, 
where J, Q, T are matrices in E?^ defined by 

I{{(l,p)) = {p,-(l), Q{{q,p)) = {Vq,p), Tip = {il),gn)ign- 

Note that 

(^,V')2 = (QV',V'')i- 

Let Al be the adjoint operator to A with respect to the corresponding scalar 
products (, k = 1,2. ft is easy to see that 

Al= -QI- aT, 

A*2 = Q-^AIQ = -IQ - aT. 

For fc = 1,2 put 

Mk = {{{Aiy9nh=o,i,...)- 

We will prove that the orthogonal complements (with respect to the correspond- 
ing scalar products) = Lq, k = 1,2. Let V G M^. Since 

+00 j 

(e*^V, 9n)k = e'^hn)k = E ^ (^*ky9n)k = (6) 

for any /; > 0, we have ip E Lq. Vice- versa, let ip G Lq. Then for any j = 0,1, . . . 
we have a chain of equalities 

= ^ (e*^V',5n)fc|t=0 = {A^llj,gn)k = (V, (^fcF^n)/^- 

It follows that tp G M^, and thus = Lq, which is equivalent to 
Let us use induction in m to prove that 

({(-42)^9nb=0,l,...,m) = {{{IQy9n}3=0.1,....rn)- 

For m = the statement is evident. The inductive hypothesis is: for some cj, 

m 

{A*2rgn = J2c,iIQygn. 

j=0 
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Applying to both sides of this equality the operator A2 = —IQ — ctF, we get 

771 

(A;)"+ig„=-^c,(/Q)-'+i5„-a((A^)™.9n,g«)2.9„e({(/Q)^5nb=o,i...,,n+i). 

Similarly, one can prove the inverse inclusion for the corresponding linear spans. 
One can prove similarly that 

L^'^ = ({(g/F5nb=o,i,...). 

It is easy to check that the following identities hold for j = 0, 1, . . . 

(/Q)2^ff„ = (-1)^(0, F^e„), 

Thus, =lyQ ly. Similarly for j = 0, 1, . . . one can check the formulas 

{Qlf'9n= (-1)^(0, ^^"e„), 

{Qlf=+'g„ = (-l)^(F^+ie„,0). 

It follows that L^^^ = y{W) ffi where V{lv) is the image of Iv under the 
mapping V . Since Iv is invariant with respect to V and V is invertible, V{lv) = 

Iv Thus, =Zy©Zy = 

Invariance with respect to the operator A follows from relations 

A(T/^e„,0r=(0,-F^"+ie„)eL(2), 
A{Q,V'e,,f = (T/^e„,0) -a(F^e„,e„)ig„ G L^^). 

The lemma is proved. □ 

Now we will use the classical result about ordinary differential equations. 
Namely, consider the following system 

2/ = /(y), yeR". (7) 

If a function v{y) G is given, we will denote its derivatives along the trajec- 
tories as 

Theorem 3.1 (Barbashin — Krasovskij's theorem [T], p. 19, Thm. 3.2) 

Let f be continuous, /(O) = and for \y\ ^ p {p > 0) there exists a (Lyapounov) 
function v{y) G such that v(0) = 0, v{y) > for y ^ 0, dv/dt ^ 0. Assume 
moreover that the set of all y for which dv/dt = does not contain any whole 
trajectory, except y = 0. Then the zero solution of ^ is asymptotically stable. 
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Lemma 3.2 The following equality holds 

Let us prove that c L^. We use Lemma [3.11 and choose v{tp) = H{i/j)- 
Then 

dv 

■^(■0) = -{Dp,p)i -Q!(?A,g„)2 ■ 

The trajectories which belong to the set of zeros of dv/dt, by definition belong 
to Lq. As Lo n L(2) = {0}^ the zero trajectory is asymptotically stable on L' '. 
This means that for any ip e i*^^^ we have H{e^'^'ip) — as i — +oo. Let us 
prove the inverse inclusion. Note that for any ip G L we have H{ip) — (ip, tp)2- 
Let € L_. We have 

Then i?(e*^V) = H{iP°) + H{iP^^^) ^ H{iP°) = by definition of It follows 
that V-''*''' = 0. Lemma [X^ and thus Theorem I2.1l are proved. 

3.2 Proof of Theorem [272] 

First, we need the following general assertion. 

Lemma 3.3 Let the spectrum ofV be simple and let {vi, . . . ,vn} be the eigen- 
vectors, which form the basis of the space . Then 

diniLo-2#{fce {!,..., iV}: Vfc ± e„}, 

where -L is in the sense o/ (, )i. 

Proof. As the subspace ly is invariant with respect to V, one can enumerate 
the vectors {vi, . . . , ujv} so that for some m 

Iv = {vi,...,Vm)- 

Thus, e„ e (wi, . . . ,w„i) = {vjrt+i, . . . jWat)^, -L in the sense of (, )i. □ 

We will calculate the spectrum of the matrix V, which acts on the vectors q 
as follows: 

{ujoqi - uji{qi^i + qi+i - 2qi), 1 < i < N - I, 
wogi - wi(g2 - gi), i = 1, 

uJoqN - i^i{qN-i - Qn), i^N. 

Lemma 3.4 The matrix V has N different eigenvalues Ao, . . . , Aat-i , given by 
Afc = Wo + 2wifl - cos^y fc 0, . . . ,iV - 1. 



N . 

The corresponding eigenvectors y''^\ . . . , have the following coordinates: 

.(k) __ Jj-l/2)k 
N 



yf^=cos7r ^-^ J ' , j = l,...,7V, fc = 0,...,iV-L 
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Lemma 13.41 can be proved by direct substitution. 
Using Lemmas 13.31 and 13.41 we get 



where 



dimio = 2#{fc e {0, . . . , TV - 1} : (y^'^) , e„) = 0} 

^2#{fce{0,...,iV-l}: cos(^M!!_M 
= 2#{fce {0,...,iV-l} : a„(fc) eZ}, 



= = 



1 / 2n - 1 



iV 



1 . 



Denoted — gcd(A^, 2n— 1). For some integers mi, rn2 such that gcd(TOi,m2) — 1, 
we have 

2n^l = dmi, N — dm2, an(k) — —( — ^-fc — 1 

2 Vm2 ^ 

As the numbers d and mi are odd, the number a^^k) is an integer if and only 
if fc = {21 — l)m2, I — 1,2, . . . From the condition fc ^ — 1 it fohows that 



I s$ - 



1 /N -1 



1 



2 V m2 

Thus 

dimio = 2 

Since {d/2N) ^ (1/2) and d is odd, we have 

d d 1 d-l 
2^ 2N 



1,7V-1 







d ■ 


= 2 








.2 ~ 


2N. 



This gives the first statement of Theorem l2.2l Statements 2 and 3 of the theorem 
follow from the results of [S], see theorems 2.1, 3.2, 4.2 therein. 
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